The PRC2.1 Subcomplex Opposes G1 Progression through Regulation of CCND1 and CCND2.

Progression through the G1 phase of the cell cycle is the most highly regulated step in cellular division. We employed a chemogenomics approach to discover novel cellular networks that regulate cell cycle progression. This approach uncovered functional clusters of genes that altered sensitivity of cells to inhibitors of the G1/S transition. Mutation of components of the Polycomb Repressor Complex 2 rescued growth inhibition caused by the CDK4/6 inhibitor palbociclib, but not to inhibitors of S phase or mitosis. In addition to its core catalytic subunits, mutation of the PRC2.1 accessory protein MTF2, but not the PRC2.2 protein JARID2, rendered cells resistant to palbociclib treatment. We found that PRC2.1 (MTF2), but not PRC2.2 (JARID2), was critical for promoting H3K27me3 deposition at CpG islands genome-wide and in promoters. This included the CpG islands in the promoter of the CDK4/6 cyclins CCND1 and CCND2, and loss of MTF2 lead to upregulation of both CCND1 and CCND2. Our results demonstrate a role for PRC2.1, but not PRC2.2, in promoting G1 progression.

Decomposition of dynamic transcriptomic responses during effector-triggered immunity reveals conserved responses in two distinct plant cell populations.

Rapid plant immune responses in the appropriate cells are needed for effective defense against pathogens. Although transcriptome analysis is often used to describe overall immune responses, collecting transcriptome data with sufficient resolution in both space and time is challenging. We reanalyzed public Arabidopsis time-course transcriptome data obtained after a low-dose inoculation of a Pseudomonas syringae strain expressing the effector AvrRpt2, which induces Effector-Triggered Immunity (ETI) in Arabidopsis. Double-peak time-course patterns were prevalent among thousands of upregulated genes. We implemented a multi-compartment modeling approach to decompose the double-peak pattern into two single-peak patterns for each gene. The decomposed peaks revealed an "echoing" pattern: the peak times of the first and second peaks correlated well across most upregulated genes. We demonstrated that two peaks likely represent responses of two distinct cell populations, which respond either cell-autonomously or indirectly to AvrRpt2. Thus, the peak decomposition extracted spatial information from the time-course data. The echoing pattern also indicated a conserved transcriptome response with different initiation times between two cell populations, despite different elicitor types. A gene set highly overlapping with the conserved gene set was also upregulated with similar kinetics during Pattern-Triggered Immunity (PTI). WRKY transcription factors appeared to underlie these similar but not identical transcriptome responses. Activation of a WRKY network via different entry-point WRKYs could explain the similar but not identical transcriptome responses elicited by different elicitor types. We discuss potential benefits of the properties of the WRKY activation network as an immune signaling network in the light of pressure from fast-evolving pathogens.

Proteome-scale movements and compartment connectivity during the eukaryotic cell cycle.

Cell cycle progression relies on coordinated changes in the composition and subcellular localization of the proteome. By applying two distinct convolutional neural networks on images of millions of live yeast cells, we resolved proteome-level dynamics in both concentration and localization during the cell cycle, with resolution of ∼20 subcellular localization classes. We show that a quarter of the proteome displays cell cycle periodicity, with proteins tending to be controlled either at the level of localization or concentration, but not both. Distinct levels of protein regulation are preferentially utilized for different aspects of the cell cycle, with changes in protein concentration being mostly involved in cell cycle control and changes in protein localization in the biophysical implementation of the cell cycle program. We present a resource for exploring global proteome dynamics during the cell cycle, which will aid in understanding a fundamental biological process at a systems level.

BridGE: a pathway-based analysis tool for detecting genetic interactions from GWAS.

Genetic interactions have the potential to modulate phenotypes, including human disease. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions; however, traditional methods for identifying them, which tend to focus on testing individual variant pairs, lack statistical power. In this protocol, we describe a novel computational approach, called Bridging Gene sets with Epistasis (BridGE), for discovering genetic interactions between biological pathways from GWAS data. We present a Python-based implementation of BridGE along with instructions for its application to a typical human GWAS cohort. The major stages include initial data processing and quality control, construction of a variant-level genetic interaction network, measurement of pathway-level genetic interactions, evaluation of statistical significance using sample permutations and generation of results in a standardized output format. The BridGE software pipeline includes options for running the analysis on multiple cores and multiple nodes for users who have access to computing clusters or a cloud computing environment. In a cluster computing environment with 10 nodes and 100 GB of memory per node, the method can be run in less than 24 h for typical human GWAS cohorts. Using BridGE requires knowledge of running Python programs and basic shell script programming experience.

RNF4 prevents genomic instability caused by chronic DNA under-replication.

Eukaryotic genome stability is maintained by a complex and diverse set of molecular processes. One class of enzymes that promotes proper DNA repair, replication and cell cycle progression comprises small ubiquitin-like modifier (SUMO)-targeted E3 ligases, or STUbLs. Previously, we reported a role for the budding yeast STUbL synthetically lethal with sgs1 (Slx) 5/8 in preventing G2/M-phase arrest in a minichromosome maintenance protein 10 (Mcm10)-deficient model of replication stress. Here, we extend these studies to human cells, examining the requirement for the human STUbL RING finger protein 4 (RNF4) in MCM10 mutant cancer cells. We find that MCM10 and RNF4 independently promote origin firing but regulate DNA synthesis epistatically and, unlike in yeast, the negative genetic interaction between RNF4 and MCM10 causes cells to accumulate in G1-phase. When MCM10 is deficient, RNF4 prevents excessive DNA under-replication at hard-to-replicate regions that results in large DNA copy number alterations and severely reduced viability. Overall, our findings highlight that STUbLs participate in species-specific mechanisms to maintain genome stability, and that human RNF4 is required for origin activation in the presence of chronic replication stress.

A scalable platform for efficient CRISPR-Cas9 chemical-genetic screens of DNA damage-inducing compounds.

Current approaches to define chemical-genetic interactions (CGIs) in human cell lines are resource-intensive. We designed a scalable chemical-genetic screening platform by generating a DNA damage response (DDR)-focused custom sgRNA library targeting 1011 genes with 3033 sgRNAs. We performed five proof-of-principle compound screens and found that the compounds' known modes-of-action (MoA) were enriched among the compounds' CGIs. These scalable screens recapitulated expected CGIs at a comparable signal-to-noise ratio (SNR) relative to genome-wide screens. Furthermore, time-resolved CGIs, captured by sequencing screens at various time points, suggested an unexpected, late interstrand-crosslinking (ICL) repair pathway response to camptothecin-induced DNA damage. Our approach can facilitate screening compounds at scale with 20-fold fewer resources than commonly used genome-wide libraries and produce biologically informative CGI profiles.

N-terminal acetylation shields proteins from degradation and promotes age-dependent motility and longevity.

Most eukaryotic proteins are N-terminally acetylated, but the functional impact on a global scale has remained obscure. Using genome-wide CRISPR knockout screens in human cells, we reveal a strong genetic dependency between a major N-terminal acetyltransferase and specific ubiquitin ligases. Biochemical analyses uncover that both the ubiquitin ligase complex UBR4-KCMF1 and the acetyltransferase NatC recognize proteins bearing an unacetylated N-terminal methionine followed by a hydrophobic residue. NatC KO-induced protein degradation and phenotypes are reversed by UBR knockdown, demonstrating the central cellular role of this interplay. We reveal that loss of Drosophila NatC is associated with male sterility, reduced longevity, and age-dependent loss of motility due to developmental muscle defects. Remarkably, muscle-specific overexpression of UbcE2M, one of the proteins targeted for NatC KO-mediated degradation, suppresses defects of NatC deletion. In conclusion, NatC-mediated N-terminal acetylation acts as a protective mechanism against protein degradation, which is relevant for increased longevity and motility.

Dimensionality reduction methods for extracting functional networks from large-scale CRISPR screens.

CRISPR-Cas9 screens facilitate the discovery of gene functional relationships and phenotype-specific dependencies. The Cancer Dependency Map (DepMap) is the largest compendium of whole-genome CRISPR screens aimed at identifying cancer-specific genetic dependencies across human cell lines. A mitochondria-associated bias has been previously reported to mask signals for genes involved in other functions, and thus, methods for normalizing this dominant signal to improve co-essentiality networks are of interest. In this study, we explore three unsupervised dimensionality reduction methods-autoencoders, robust, and classical principal component analyses (PCA)-for normalizing the DepMap to improve functional networks extracted from these data. We propose a novel "onion" normalization technique to combine several normalized data layers into a single network. Benchmarking analyses reveal that robust PCA combined with onion normalization outperforms existing methods for normalizing the DepMap. Our work demonstrates the value of removing low-dimensional signals from the DepMap before constructing functional gene networks and provides generalizable dimensionality reduction-based normalization tools.

RNF4 and USP7 cooperate in ubiquitin-regulated steps of DNA replication.

DNA replication requires precise regulation achieved through post-translational modifications, including ubiquitination and SUMOylation. These modifications are linked by the SUMO-targeted E3 ubiquitin ligases (STUbLs). Ring finger protein 4 (RNF4), one of only two mammalian STUbLs, participates in double-strand break repair and resolving DNA-protein cross-links. However, its role in DNA replication has been poorly understood. Using CRISPR/Cas9 genetic screens, we discovered an unexpected dependency of RNF4 mutants on ubiquitin specific peptidase 7 (USP7) for survival in TP53-null retinal pigment epithelial cells. TP53-/-/RNF4-/-/USP7-/- triple knockout (TKO) cells displayed defects in DNA replication that cause genomic instability. These defects were exacerbated by the proteasome inhibitor bortezomib, which limited the nuclear ubiquitin pool. A shortage of free ubiquitin suppressed the ataxia telangiectasia and Rad3-related (ATR)-mediated checkpoint response, leading to increased cell death. In conclusion, RNF4 and USP7 work cooperatively to sustain a functional level of nuclear ubiquitin to maintain the integrity of the genome.

CD69 marks a subpopulation of acute myeloid leukemia with enhanced colony forming capacity and a unique signaling activation state.

In acute myeloid leukemia (AML), leukemia stem cells (LSCs) have self-renewal potential and are responsible for relapse. We previously showed that, in Mll-AF9/NRASG12V murine AML, CD69 expression marks an LSC-enriched subpopulation with enhanced in vivo self-renewal capacity. Here, we used CyTOF to define activated signaling pathways in LSC subpopulations in Mll-AF9/NRASG12V AML. Furthermore, we compared the signaling activation states of CD69High and CD36High subsets of primary human AML. The human CD69High subset expresses low levels of Ki67 and high levels of NFκB and pMAPKAPKII. Additionally, the human CD69High AML subset also has enhanced colony-forming capacity. We applied Bayesian network modeling to compare the global signaling network within the human AML subsets. We find that distinct signaling states, distinguished by NFκB and pMAPKAPKII levels, correlate with divergent functional subsets, defined by CD69 and CD36 expression, in human AML. Targeting NFκB with proteasome inhibition diminished colony formation.

Immunophenotypically-defined murine AML stem cells harbor self-renewing and non-self-renewing subsets that display unique signaling characteristics.CD69, an NFκB target gene, marks a subset of human AML with increased colony forming capacity and reduced proliferation.NFκB activation correlates with the global signaling pathway activation state in human AML.

Reproducibility metrics for context-specific CRISPR screens.

CRISPR screens are used extensively to systematically interrogate the phenotype-to-genotype problem. In contrast to early CRISPR screens, which defined core cell fitness genes, most current efforts now aim to identify context-specific phenotypes that differentiate a cell line, genetic background, or condition of interest, such as a drug treatment. While CRISPR-related technologies have shown great promise and a fast pace of innovation, a better understanding of standards and methods for quality assessment of CRISPR screen results is crucial to guide technology development and application. Specifically, many commonly used metrics for quantifying screen quality do not accurately measure the reproducibility of context-specific hits. We highlight the importance of reporting reproducibility statistics that directly relate to the purpose of the screen and suggest the use of metrics that are sensitive to context-specific signal. A record of this paper's transparent peer review process is included in the supplemental information.

Fanconi anemia-associated chromosomal radial formation is dependent on POLθ-mediated alternative end joining.

Activation of the Fanconi anemia (FA) pathway after treatment with mitomycin C (MMC) is essential for preventing chromosome translocations termed "radials." When replication forks stall at MMC-induced interstrand crosslinks (ICLs), the FA pathway is activated to orchestrate ICL unhooking and repair of the DNA break intermediates. However, in FA-deficient cells, how ICL-associated breaks are resolved in a manner that leads to radials is unclear. Here, we demonstrate that MMC-induced radials are dependent on DNA polymerase theta (POLθ)-mediated alternative end joining (A-EJ). Specifically, we show that radials observed in FANCD2-/- cells are dependent on POLθ and DNA ligase III and occur independently of classical non-homologous end joining. Furthermore, treatment of FANCD2-/- cells with POLθ inhibitors abolishes radials and leads to the accumulation of breaks co-localizing with common fragile sites. Uniformly, these observations implicate A-EJ in radial formation and provide mechanistic insights into the treatment of FA pathway-deficient cancers with POLθ inhibitors.

Systems-level transcriptional regulation of Caenorhabditis elegans metabolism.

Metabolism is controlled to ensure organismal development and homeostasis. Several mechanisms regulate metabolism, including allosteric control and transcriptional regulation of metabolic enzymes and transporters. So far, metabolism regulation has mostly been described for individual genes and pathways, and the extent of transcriptional regulation of the entire metabolic network remains largely unknown. Here, we find that three-quarters of all metabolic genes are transcriptionally regulated in the nematode Caenorhabditis elegans. We find that many annotated metabolic pathways are coexpressed, and we use gene expression data and the iCEL1314 metabolic network model to define coregulated subpathways in an unbiased manner. Using a large gene expression compendium, we determine the conditions where subpathways exhibit strong coexpression. Finally, we develop "WormClust," a web application that enables a gene-by-gene query of genes to view their association with metabolic (sub)-pathways. Overall, this study sheds light on the ubiquity of transcriptional regulation of metabolism and provides a blueprint for similar studies in other organisms, including humans.

Dimensionality reduction methods for extracting functional networks from large-scale CRISPR screens.

CRISPR-Cas9 screens facilitate the discovery of gene functional relationships and phenotype-specific dependencies. The Cancer Dependency Map (DepMap) is the largest compendium of whole-genome CRISPR screens aimed at identifying cancer-specific genetic dependencies across human cell lines. A mitochondria-associated bias has been previously reported to mask signals for genes involved in other functions, and thus, methods for normalizing this dominant signal to improve co-essentiality networks are of interest. In this study, we explore three unsupervised dimensionality reduction methods - autoencoders, robust, and classical principal component analyses (PCA) - for normalizing the DepMap to improve functional networks extracted from these data. We propose a novel "onion" normalization technique to combine several normalized data layers into a single network. Benchmarking analyses reveal that robust PCA combined with onion normalization outperforms existing methods for normalizing the DepMap. Our work demonstrates the value of removing low-dimensional signals from the DepMap before constructing functional gene networks and provides generalizable dimensionality reduction-based normalization tools.

Genetic architecture of protein and oil content in soybean seed and meal.

Soybean is grown primarily for the protein and oil extracted from its seed and its value is influenced by these components. The objective of this study was to map marker-trait associations (MTAs) for the concentration of seed protein, oil, and meal protein using the soybean nested association mapping (SoyNAM) population. The composition traits were evaluated on seed harvested from over 5000 inbred lines of the SoyNAM population grown in 10 field locations across 3 years. Estimated heritabilities were at least 0.85 for all three traits. The genotyping of lines with single nucleotide polymorphism markers resulted in the identification of 107 MTAs for the three traits. When MTAs for the three traits that mapped within 5 cM intervals were binned together, the MTAs were mapped to 64 intervals on 19 of the 20 soybean chromosomes. The majority of the MTA effects were small and of the 107 MTAs, 37 were for protein content, 39 for meal protein, and 31 for oil content. For cases where a protein and oil MTAs mapped to the same interval, most (94%) significant effects were opposite for the two traits, consistent with the negative correlation between these traits. A coexpression analysis identified candidate genes linked to MTAs and 18 candidate genes were identified. The large number of small effect MTAs for the composition traits suggest that genomic prediction would be more effective in improving these traits than marker-assisted selection.

Survival-based CRISPR genetic screens across a panel of permissive cell lines identify common and cell-specific SARS-CoV-2 host factors.

SARS-CoV-2 depends on host cell components for infection and replication. Identification of virus-host dependencies offers an effective way to elucidate mechanisms involved in viral infection and replication. If druggable, host factor dependencies may present an attractive strategy for anti-viral therapy. In this study, we performed genome wide CRISPR knockout screens in Vero E6 cells and four human cell lines including Calu-3, UM-UC-4, HEK-293 and HuH-7 to identify genetic regulators of SARS-CoV-2 infection. Our findings identified only ACE2, the cognate SARS-CoV-2 entry receptor, as a common host dependency factor across all cell lines, while other host genes identified were largely cell line specific, including known factors TMPRSS2 and CTSL. Several of the discovered host-dependency factors converged on pathways involved in cell signalling, immune-related pathways, and chromatin modification. Notably, the chromatin modifier gene KMT2C in Calu-3 cells had the strongest impact in preventing SARS-CoV-2 infection when perturbed.

Scalable CRISPR-Cas9 chemical genetic screens in non-transformed human cells.

Pooled lentiviral CRISPR-Cas9 screens are utilized for assessing the differential sensitivity or resistance of many single-gene knockouts to a compound. Here, we present a scalable approach for high-throughput compound screening by utilizing a small custom library. We describe steps to perform a proof-of-principle chemical screen in non-transformed hTERT RPE-1 TP53-/- cells with higher coverage and greater timepoint resolution compared to genome-wide screens. This approach can be adapted for use in various cell lines, compounds, and other focused sgRNA libraries.

Repression of essential cell cycle genes increases cellular fitness.

A network of transcription factors (TFs) coordinates transcription with cell cycle events in eukaryotes. Most TFs in the network are phosphorylated by cyclin-dependent kinase (CDK), which limits their activities during the cell cycle. Here, we investigate the physiological consequences of disrupting CDK regulation of the paralogous repressors Yhp1 and Yox1 in yeast. Blocking Yhp1/Yox1 phosphorylation increases their levels and decreases expression of essential cell cycle regulatory genes which, unexpectedly, increases cellular fitness in optimal growth conditions. Using synthetic genetic interaction screens, we find that Yhp1/Yox1 mutations improve the fitness of mutants with mitotic defects, including condensin mutants. Blocking Yhp1/Yox1 phosphorylation simultaneously accelerates the G1/S transition and delays mitotic exit, without decreasing proliferation rate. This mitotic delay partially reverses the chromosome segregation defect of condensin mutants, potentially explaining their increased fitness when combined with Yhp1/Yox1 phosphomutants. These findings reveal how altering expression of cell cycle genes leads to a redistribution of cell cycle timing and confers a fitness advantage to cells.

Context-dependent regulation of ferroptosis sensitivity.

Ferroptosis is an important mediator of pathophysiological cell death and an emerging target for cancer therapy. Whether ferroptosis sensitivity is governed by a single regulatory mechanism is unclear. Here, based on the integration of 24 published chemical genetic screens combined with targeted follow-up experimentation, we find that the genetic regulation of ferroptosis sensitivity is highly variable and context-dependent. For example, the lipid metabolic gene acyl-coenzyme A (CoA) synthetase long chain family member 4 (ACSL4) appears far more essential for ferroptosis triggered by direct inhibition of the lipid hydroperoxidase glutathione peroxidase 4 (GPX4) than by cystine deprivation. Despite this, distinct pro-ferroptotic stimuli converge upon a common lethal effector mechanism: accumulation of lipid peroxides at the plasma membrane. These results indicate that distinct genetic mechanisms regulate ferroptosis sensitivity, with implications for the initiation and analysis of this process in vivo.